DYNAMICS OF SOLITONS AND QUASISOLITONS OF CUBIC THIRD-ORDER 

NONLINEAR SCHRODINGER EQUATION 
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The dynamics of soliton and quasisoliton solutions of cubic third order nonlinear Schrodinger 
equation is studied. The regular solitons exist due to a balance between the nonlinear terms and 
(linear) third order dispersion; they are not important at small is the coefficient in the 

third derivative term) and vanish at as 0. The most essential, at small as, is a quasisoliton 
emitting resonant radiation (resonantly radiating soliton). Its relationship with the other (steady) 
quasisoliton, called embedded soliton, is studied analytically and in numerical experiments. It is 
demonstrated that the resonantly radiating solitons emerge in the course of nonlinear evolution, 
which shows their physical significance. 



PACS numbers: 42.65.Tg, 52.35. Mw 



I. INTRODUCTION 

Equations describing soliton processes are usually ob- 
tained by certain approximating procedures affecting 
nonlinearity and dispersion. Now, in extensive studies 
of ultrafast processes, the classical approximations often 
appear to be insufficient and higher order effects become 
of importance. A typical example is high speed systems 
like nonlinear transmission lines in femtosecond regime 
for soliton communications, etc. For such systems soliton 
solutions in classical sense do not generally exist. Only 
for very specific choices of parameters one can find lo- 
calized solutions. Instead of regular solitons there may 
appear nonlocal steady or (and) unsteady soliton like 
structures (which may be called quasisolitons). These 
may, however, have significant importance for the non- 
linear dynamics. 

In this paper we consider dynamics described by the 
extended third-order cubic nonlinear Schrodinger (NLS) 
equation 



1 



|$|2^ 



+1038x^1 



0, 



(1) 



which plays an important role in many nonlinear prob- 
lems, in partiqular, in the nonlinear fiber optics [1,2]. The 
third derivative term describes higher order correction to 
the linear dispersion and the cubic terms with the first 
derivative are dispersive nonlinear corrections. Here, we 
assume that all a„ ( 1 < ctji < 3) are real parameters. 

Eq.(l) has regular soliton solutions [3-5], vanishing at 
|X| 0. In particular cases, Eq.(l) may have exact A^- 
soliton solutions [4] or even be integrable [5]. An impor- 
tant feature of the solutions describing regular solitons is 
that they are degenerating at as — s- (see below), which 
substantially reduces their physical importance. Their 
numerical investigation at ai^2 1 where done e.g. in 
Ref.[6] (see also references therein) 



There also exist some solutions to Eq.(l), describing 
quasisolitons. One of them is a steady solution looking 
like a soliton type pulse embedded into a small ampli- 
tude plane wave (a soliton on plane wave pedestal). We 
call it embedded soliton (ES), using the word proposed in 
a different context in Ref. [7]. Apart from that, Eq. (1) 
has other type of quasisoliton solutions which describe 
soliton-like pulses permanently emitting resonantly gen- 
erated radiation (see, e.g., Refs. [8-12] for the particular 
case at ai = a2 = and [13,14] for the fuU Eq.(l)). We 
will call them resonantly radiating solitons (RRS); they 
are unsteady because of losses caused by the radiation. 
The lifetime of radiating soliton is sufficiently large if 03 
is small enough and, naturally, one can speak about the 
soliton only in this quasisteady case. (However, at long 
times, the losses caused by the radiation become essential 
for applications [2].) At a3 ^ 0, the RRS turns into the 
regular soliton of Eq.(l) without the third derivative. 

There is an interesting and important connection be- 
tween the two types of quasisolitons, ES and RRS, and 
we demonstrate it in numerical simulations. We will also 
study the role of solitons and quasisolitons in the non- 
linear processes described by Eq.(l). In partiqular, the 
regular solitons and RRS compete between themselves 
in nonlinear processes. As far as the regular solitons dis- 
appear at as ^ and the RRS have short lifetimes at 
large as, it is clear that at large 0:3 the regular solitons 
are more important while the RRS may play a decisive 
role at sufficiently small a^. The most interesting case 
is, of course, small 0:3 , because the third derivative term 
emerges as the result of an expansion. (The effect of next, 
4th derivative, term can be seen in Refs. [13,14]; at cer- 
tain relalionship between the coefRcients before third and 
fourth derivatives there can be no radiation at all which, 
naturally, may happen also when other high-order terms 
are taken into account.) 

The paper is organized as follows. In Sec. I we de- 
scribe important properties of Eq.(l), which are used 
below. In particular, we discuss regular soliton solu- 
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tions [3] , Galilean transformations and conservation laws 
[13,14]. The embedded solitons are studied in Sec. 3. As 
we have already mentioned, the ES consists of a steady 
soliton-like pulse on the plane wave background. This 
structure is rater common in different nonlinear highly 
dispersive systems [9-20]. At as — > , the plane wave 
disappears and the pulse turns into a regular soliton of 
Eq. (1) at as = . The plane wave amplitude, in- 
creasing with 0:3 , may become ^instable at sufficiently 
large as due to the modulational instability. Consider- 
ing small as, we show that the pulse part of ES is rather 
close to the pulse in RRS and the wavenumber of the 
plane wave coincides with the wavenumber of resonantly 
emitted radiation (by the RRS) . Next, we introduce the 
cut off operation, cutting off both wave wings from the 
ES (Sec. IV). Then we see that the remained soliton-like 
pulse is transformed into the RRS, emitting radiation 
only in one direction , according to the direction of the 
group velocity of resonant radiation. The resonant radi- 
ation disappears at ai = 60:3. In this case, considered 
in Sec.V, we show that Eq.(l) can be transformed, by 
means of the Galilean transformation, to the complex 
modified Korteweg-de Vries (MKdV) equation . At par- 
ticular initial conditions, it is reduced to the real MKdV 
equation which is integrable and therefore has N- soliton 
solutions [22] . This does not mean the complete integra- 
bility of Eq.(l) at ai = Sas , because the reduction to 
the real MKdV equation is possible only for particular 
initial conditions. However, the Painleve analysis [21] 
of the above mentioned complex MKdV equation shows 
that at 60:3 = ai = 2a2 it possesses the Painleve prop- 
erty. This is just the case when Eq. (1) is integrable 
[5]. If a2 = 0, we arrive at the case considered by Hirota 
[4] who has shown the existence of - soliton solutions 
in this case. Our approach, however, is more general 
and (in fact, as in Ref. [23] ) we can see that the N 
-soliton solutions are possible even for nonvanishing a2 
(at a2 — 3a3 ). Then we present a numerical study 
of the initial value problem in the Hirota casc:ai = Gas, 
a2 = taking as initial conditions complex pulses and 
one can see how the robust Hirota solitons are emerg- 
ing from these pulses. In Sec. VI, we describe numerical 
experiments at ai ^ 6as. Solving the initial value prob- 
lem, we show that, at sufficiently small as, the initial 
disturbances decay into the resonantly radiating solitons 
(RRS). On the other hand, at small as, the Potasek- Ta- 
bor solitons [3] were not detected. This indicates that the 
RRS are important physical objects playing a significant 
role in nonlinear dynamics. In Sec. VII , a summary of 
obtained results is given. In Appendix A, the modula- 
tional instability of a plane wave, described by Eq.(l), is 
considered; it is helpful in the study of the stability of 
embedded solitons. In Appendix B we investigate ana- 
lytically, by means of conservation laws, the evolution of 
RRS caused by the radiation. This analysis is in agree- 
ment with numerical results described in Sees. IV and VI. 

II. IMPORTANT PROPERTIES OF EQ. (1) 



a. Exact soliton solutions 

First, we discuss the exact soliton solutions of Eq.(l). 
They can be written in the form 

= o sech[b{X - V,T)] e^''^'^'^ , (2a) 

^2^_6a3 ^2 (c,,^_2a2), (26) 

ai -|- 2a2 

Vs = K — 3asK^ -|- as6^, (2c) 

uj= -K^ - asK'' 6^ {2d) 

ai + 2a2 - 6as , 

K= . (2e) 

12a2a3 ^ ^ 

In fact, these are the soliton solutions found by Potasek 
and Tabor [3] (with corrected misprints). We will call 
Eqs. (2a-e) Potasek- Tabor (PT) soliton solutions. At 
as = 0, PT solitons do not exist: this signifies that their 
existence is a result of a ballance between the nonlin- 
ear terms and the linear third order dispersive term . If 
a2 ^ 0, from Eq. (2e) it follows that k is a fixed num- 
ber, uniquely determined by the coefficients of Eq.(l). 
However, the solitons (2) exist even at a2 = provided 
that 

ai = Gas. (3) 

In this case k can be arbitrary , because Eq. (2e) follows 
from the equation 12a2asK = ai — 6as -|- 2a2. 

In case (3) and a2 = 0, Eq. (1) is the so called Hirota 
equation that can be transformed to the complex mod- 
ified Korteweg-de Vries equation which has N -soliton 
solutions [4]. For the Hirota soliton, from Eq.(2b) it fol- 
lows 

b = a. (4) 

In the other particular case, 603 = ai = 2a2, Eq. (1) is 
integrable [5]. (See also Sec.V.) Some results for solitons 
(2) with a2 ^ were reported in Ref. [6] (see also refer- 
ences therein) . In particular, it was shown numerically 
that they emerge in a solution of initiall value problem. 
However, this result was obtained at ai = as = 1; at 
small as, as it will be shown below, they cannot compete 
with RRS. 

Apart from the solution (2a-e), describing "bright" 
solitons, which was called in Ref. [3] sec/i- family, Potasek 
and Tabor have also found tanh-iamily which will not be 
discussed here. (Some generalizations, with both sech 
and tanh terms, are derived in Ref. [27] by means 
of rather tedious algebra; the results of this paper are 
mostly contained in Ref. [3] and the generalizations [27], 
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which arc vaHd for very specific choices of the parame- 
ters, can be obtained by the much simpler approach of 
Potasek and Tabor.) The results of Ref.[3] are repeated 
also in some other papers (e.g. [23]). 

b. Galilean transformation 

In the following we shall frequently use the Galilean 

transformation describing the transition to the reference 
frame moving with the velocity V . For Eq. (1) it reads 



^{X, T) = ipiX - VT, T) ex.p[i{KX - QT)] 
where K and Q. are defined by equations 
V = K- 3a3i^'^ 

2 

The function ip{x,t) with 

x = X-VT , t = T 
satisfies the equation [13,14] 
1 



(5) 

(6) 
(7) 

(8) 



with 



+ia3dlip = 



a2 = 1 - 6a3K, 



1 - aiK. 



(9) 



(10) 



(11) 



c. Conservation laws 

One can check by means of straightforward calcula- 
tions that Eq. (1) has the following integral of motion 



N 



/OO 
\^{X, 
-OO 



T)fdX . 



At 



Ct2 = 0, 

we have two other conserving integrals [24, 13,14] 



(*ax** - ^*dx'^)dX , 



(12) 



(13) 



(14) 



/OO 1 1 1 

{-\dx^\^--\^\^--iai\^\\^*dx^-^dx^* 



-m3(**9^* - ^d^x^*)}dX. 



(15) 



III. EMBEDDED SOLITONS 

We start with the solutions to Eq. (1) of the form 

vI/(X,r) = x(x)e'^*, (16) 

where x and t are defined in (8). Substituting (16) into 
Eq.(l), we arrive at the ordinary differential equation 



-iVd^X+ 2^lx+\x\'^X + i<^i\x\'^da;X + 'ia2Xdx\x\'' 



(17) 



Imposing periodic boundary conditions at the ends of a 
sufficiently broad interval and considering Fand x(0) as 
given parameters and A as the eigenvalue, we arrive at a 
nonlinear eigenvalue problem which can be numerically 
solved by a kind of shooting method. An example of such 
solution at a2 = and 

ai = -0.3 , as = -0.1 , x(0) = 1.5 , V = (18) 

is shown in Fig.l(a,b) . It is a nonlocal steady pulse with 
I ^ (a;) I max = |^(a;o)| and small symmetric "wings". At 
large |a; — xoj we can linearize Eq. (17); then we see that 

X{x) « const e'-^^ (ja; - | » 1) , (19) 

where k is a real root of cubic algebraic equation 

A- kF + (1/2)^2 _ q,3k3 = . (20) 

Finding A from the numerical solution of the eigenvalue 
poblem and k from (20), we have for the case (18) 



A = 0.979 , 



-5.34 



(21) 



On the other hand we can determine k directly from x{^) 
which is found numerically from Eq. (17), together with 
A. From Fig. lb we see that the numerically found 
X{x) has indeed asymptotic behavior (19) with k « —5.3 
which is very close to the root of Eq. (20), written in 
(21). This agreement is an evidence of the correctness 
of the numerical solution of eigenvalue problem by the 
shooting method. In a similar way, the solution of eigen- 
value problem (17) and Eq. (20) for a2 = and 

ai = -0.8 , as = -0.1 , x(0) = 1-5 , V = -0.225 

(22a) 

gives 



A= 1.209 



-5.03 . 



(226) 
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FIG. 1. Numerical solution of Eq. (17) at parameters 
(18). (a) X versus x; solid line: IxC^^)! = l*(*)l; dashed 
line: Rex{x). (b) argx versus x; at sufficiently large |a; — a;o|, 
argxix) is a linear function: argx ^ 'A-'^ which permits 

to measure k (k ^ —5.2). 

The asymptotic behavior of numerically found x{x) 
gives the same «; ~ —5.03. 

For the other case with a2 = and 

ai = -1.1 , Q!3 = -0.1 , x(0) = 1.5 , V = 0.35 (23a) 

we have 

A = 1.937, K=-6.12. (236) 

while from the numerically found x{x) we obtain n = 
-6.1. 

In fact, similar quasisoliton solutions, looking like 
embedded solitons, were obtained, by different ap- 
proaches, for other highly dispersive systems as well [15- 
17,9,11,12,17-20]. 

Now consider the embedded solitons, starting from 
Galilean transformation (5). Writing it in the form 

*(X, T) = il,{x, t) e'^'-" eyi-p[i{KV - ^)t] (24) 

and assuming that 

1 



(25) 



where A is a constant parameter, we compare (24) with 
(16). Then we have 



tp{x) = x{x)exp{-iKx) , 



X'^ = 2{A + n- KV) . 



(26) 



(27) 



From (25) and (9) we arrive at the following equation 
for V'(a^) 



lx2T , 1 ,2 



+iaidl'4) = , 



(28) 



Substituting here (26), (27) and taking into account (6), 
(7) and (10), (11) we have Eq.(17) as one should expect. 



From Eq. (26) it is seen that < 1 at \x\ > 1. 
Linearizing Eq. (28) we obtain 



\x\ » 1), 



(29) 



where is a root of the equation 

2a3fc^ - a2k'^ - A^ = . (30) 

Substituting (19) and (29) into Eq. (26), we see that 

K = k + K. (31) 

Then, using (6), (7) and (10), (27) we can easily prove 
that Eqs. (30) and (20) arc equivalent, as it should bo. 
At small as, the solution of Eq.(28) can be written as 



-ipix) = [us{x) + f{x)]ex.p[i(j)s{x)]. 



(32) 



with small ,f{x). Here Us{x) and (l)s{x) are defined by the 
requirement that 

F{x) = Us{x)exp[i<l)s{x)] (33) 

is a soliton solution of Eq. (28) without the last term i.e. 

~X^F+^a2dlF+q\F\''F+iai\F\^d,F+ia2Fd,\F\'' = 

(34) 

and F{x) ^ at a; = ±oo. Solving Eq. (34), we have 



Us{x) = Aw — 
V Q 



cosh ( -^^x 



-1/2 



(35) 



(j)s{x) 



«! -I- 2a2 
2A 



arctan 



tanh 



a2q 



V4A2A2 + a2q^ ' 



^2 _ 4ai(ai + 2a2) - (ai + 2a2f 
12 



(36) 

(37) 



(38) 



Eqs. (35)- (38) were obtained in Ref. [14] for a pulse part 
of RRS. From (35) it follows that the soliton amplitude 

is 



uq = Ai 



2p 



_ (i+p)g 

and its width is given by 

5 = ^faijX . 



(39) 



(40) 



Thus one must require 02 > and from Eq. (10) it 
follows that at a-^K > 0, 



1^1 < glasj 



(41) 
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which is a restriction on the soHton velocity (6). 

A small term f{x) in Eq.(32) expresses the effect of 
the third order dispersion. At \x\ y/{a2)/\, or less, it 
describes a modification of the pulse, arising due to the 
last term in Eq. (28), while at large x 



f{x) « ijjix) e 



ikx 



(42) 



where fc is a root of Eq. (30). Note, that this equation 
was derived in Ref. [14] for the wavenumber of the res- 
onantly generated radiation by the RRS in the reference 
frame where the RRS is at rest. From all that we con- 
clude that (32) coincides with the asymptotic expression 
for the RRS and its radiation at large t. 

Let us now compare the embedded solitons, obtained 
above numerically, with the solution (32) at condition 
(13) . From Eq.(38) we have 



Eqs. (6) and (10) give 



K 



1 - VI - 1203^ 
6a3 



02 = \/l - Ua^V 
Using (39) and (37), wo obtain 



Un = 



(43) 



(44) 



(45) 



(46) 



~ 2^2 ' 

where q is defined in (11). From (27) and (6), (7) we find 

1 



2 ( A - ^K^ 



(47) 



Then we have for the case (18), (21), which is shown in 
Fig. 1: uo « 1.37. For the case (22) uq « 1.41, and for 
(23) Uq « 1.42. So, in all three cases we approximately 
have Uq 1.4 which is rather close to x(G) = 1.5, as- 
sumed in the nonlinear eigenvalue problem for all three 
cases. This supports the conclusion that expressions (35) 
and (36) approximately describe the pulse in the embed- 
ded soliton (at small ). 

Now consider f{x). At small as, the roots of Eq. (30) 
have simple analytical expressions [14]. Neglecting the 
first term in (30), we have two smallest roots 



±i- 



X 



Substituting this into (42), we have 

•ip w f{x) exp(=FA/ v^a;) . 



(48) 



(49) 



This is in agreement with the asymptotic behavior of ex- 
pression (33) at large x. The third root can be approx- 
imately obtained if one neglects the last term in (30). 
This gives 



k w a2/2a3 , 



(50) 



which approximately determines the wavenumber of the 
plane wave in the wings. Expressions (50) and (31) ap- 
proximately give the roots of Eq. (20) in analytical form. 
From (50) it follows 



sgn k = sgn 0:3 . 



(51) 



And, finally, from the results of Ref. [14] it follows that 
at small as 



fix) « B 



a2\k\ 



exp 



4A 



where 5 is a complex constant with \B 
expression is valid at 



a2|A;| 



> 1 



1 H — arcsinp 

(52) 

I ~ 1 - 10. This 



(53) 



i.e. when the wavenumber k is much larger than the 
inverse width of the soliton (33). In expression (52), it 
is also assumed that 2^A A/Va2|fci,2| (this does not 
exclude AA < 1). If 

2AX ~ A/V^|fci,2| (or 2AX < A/V^|fci,2|), (54) 
which may be satisfied only at p w 1 and g w 1, we have 

TT^a^\kl 



f{x) = B./a^\k\e-xp 



2A 



Akx 



(55) 



It is easy to check that we can arrive at Eq. (55) by 
substituting in Eq. (52) the first of conditions (54) and 
p=l. 

From the all foregoing, one can see a connection be- 
tween the embedded and resonantly radiating solitons. 
In the next section we present numerical experiments dis- 
closing this connection from another side. 

IV. THE CUT OFF OPERATION 

Let us define the cut off operation transforming the 
embedded soliton into an isolated pulse. Turning to the 
function x(x) in (16), we write 



Xcut{x) = x{x)r{x) 



(56) 



where r{x) is a cutting factor which we take in the form 



r{x) = 



tanh 



X — xo + Ax 
7 



tanh 



X — Xo — Ax 

7 



(57) 

Here xo is the center of the pulse and Ao; > is the 

width of the cutted function Xcut{x). According to (57), 
r{x) vanishes at |a; — xq] 00 and the positive parameter 
7 characterizes the "sharpness" of vanishing. Assuming 
that 7 is small enough and Ax is such that the factor 
r{x) cuts off only the wings without essential disturbing 
the pulse in x{x) ( 7 and Ax can be properly chosen in 
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numerical tests) we then take the pulse Xcut{x) as the 
initial condition to Eq.(l) (see Fig. 2a) 



^{X,0) = Xcut{X) . 



(58) 



In Fig. 2b one can sec that the cutted pulse emits radia- 
tion at t > 0. However, on the left hand side of the pulse 
the radiation spreads out with time; this shows that on 
the left hand side there is a transient radiation, emitted 
at small t due to initial condition (58). On the contrary, 
on the right hand side we see a wave train with approx- 
imately constant amplitude, with the front propagating 
to the right; so, the length of the wavetrain increases with 
time. Therefore the cutted pulse permanently emits ra- 
diation to the right. The spectrum, shown in Fig.lc, has 
a peak at k approximately equal to the root of Eq. (20) 
with account of finite velocity of the pulse at t = 138 
(it is still rather close to k from Eq. (21), i.e. to the 
wavenumber of the wing waves in embedded soliton). 

Analyzing the time behavior of arg^E", we find that 
d{aig'i/) / dt — A{t) is a slow function of t, with A(0) « 
0.98 and A(130) w 0.85. Note that this A(0) coincides, 
with a good accuracy, with (21). Similar results were 
obtained for cutted pulses in cases (22) and (23). 

We conclude that the cutted pulses are radiating soli- 
tons; and the radiation is permanently emitted only in 
one direction. As far as k is connected with k from Eq. 
(31), which is the real root of Eq. (30), obtained from 
the resonant condition, we conclude that the cut off op- 
eration transforms ES ino RRS. The front of radiated 
wavetrain propagates with the group velocity U {k) given 
by [14] 



U{k) = a^k — 3a3fc^ « — a2/4a3 . 



(59) 



is in agreement with the results presented in Figs. 2b 
and 2c. 

To describe analytically the whole system after the cut 
off, one can use at small az the equation [14] 

i}{x,t) = [us{x) + ri{x,t)]ex.p[i(j)s{x)] , (61) 

where Ua{x) and (psix) are given (in adiabatic approxi- 
mation) by Eqs. (35) and (36) and //(x, y), at large x and 
t, has the following asymptotic expression 



ri{x,t)^ f{x)Q{Ux)@{\U\t-\x\) . 



(62) 



Here f{x) is given by Eq. (52) and Q{Y) is the Heaviside 
function 



e(F) = 1 (y > 0) , Q{Y) = (F < 0) 



(63) 



Eq. (62) expresses that the soliton radiates in the direc- 
tion of group velocity U and the radiation front propa- 
gates with the velocity \U\. The mentioned above adia- 
batic approximation means that at sufRciently small 0:3, 
the radiation is so small that the soliton parameters in 
Us{x) and Vs(a;) (as well as the wavenumber k) can be 
considered as constant. However, the soliton losses may 
be essential at large times. 



FIG. 3. Temporal behavior of the radiating soliton at 
Qi = —0.3, a2 = 0, as = —0.1. (a) |*|max versus T=t ; 
(b) the soliton position versus t. 



FIG. 2. Numerical solution of Eq. (1) at initial condition 
(58) and ai = —0.3, 02 = 0, 03 = —0.1. (a) The initial 
cutted pulse. Full line: I*!'! versus x&tT=0; dotted line depicts 
the corresponding embedded soliton. (b) The cutted pulse at 
T=138 ; full line: |^| versus x\ dashed line: Re^ versus x . 
{c)logiQ\^k\ versus k( '^k is the Fourier transform of '^{x); 
one can see a peak near the root of resonance equation (20). 



Prom this it follows that 

sgnU{k) = — sgnfc = — sgnas 



(60) 



Therefore at 03 < it should be A; < and U{k) > 0. 
This means that the soliton in Fig. 2a should perma- 
nently emit radiation to the right while the peak in the 
spectrum of the wavetrain should be at negative k. This 



The variation of soliton parameters, caused by the 
radiation, can be estimated by means of the integrals 
of motion (Appendix B). A vast information about the 
soliton evolution, caused by the radiation, follows from 
the munerical solution of Eq.(l). For the parameters 
used in Figs. 2, the results following from this solu- 
tion are presented in Fig. 3. One finds that the mean 
value of l^'lmax = uo{t), shown in Fig. 3a, decreases 
logarithmically, similar to what was found for the case 
ai = a2 = [25,26]. Analyzing the soliton position 
(Fig. 3b) one can find that the soliton velocity V{t) in- 
creases from V{0) = to 1/(138) « 0.29, i.e. the soliton 
is accelerating and the function V{t) increases (also log- 
arithmically). Now, using A(0) w 0.98, A(138) w 0.85 
and Eqs. (47), (44), wc have A(0) « 1.4, A(138) « 1.28. 
Therefore, the soliton width A = sJa^jX increases from 
A(0) « 0.71 to A(138) « 0.99, i.e. the radiating soliton 
is widening with time . All that can be considered as a 
numerical confirmation of analytical results obtained in 
Appenix B from the conservation laws, in particular that 



sgn-^ = — sgnfc = sgnt/ 



(64) 
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which means that the soHton is accelerating in the direc- 
tion of the group velocity of resonant radiation. 

V. SPECIAL CASE ai = 603- HIROTA SOLI- 
TONS 

The coefficient q in Eq. (9) disappears at 

K=l/ai. (65) 

Substituting this in Eq. (10) , wc have a2 = 1 — Qa^/ai. 
Therefore, at condition (3) the coefficient 02 also vanishes 
and Eq. (9) takes the form 

9tV + ^aM^d^^ + a2V'5.|^|' + azdl^ = . (66) 

which is related to Eq.(l) by a Galilean transformation 
(5), (8) with 

^=1^' ^=6^' "=lokf 
Consider a particular solution of Eq.(66) 

i>{x,t)^e''ax,t) , (68) 

where = const and 9, C{x,t) are real. Then Eq. (66) 
is reduced to the modified Korteweg-de Vries (MKdV) 
equation 

dtC + 2(3^3 + a2)C^9xC + a3dlC = • (69) 

It is completely integrable and, in particular, has exact 
N - soliton solutions [22] if 02 7^ ~3a^ = — (l/2)ai [cf. 
with Eq. (2b)]. As far as Eq. (66) is reduced to Eq. 
(69) only for particular initial conditions, it may not be 
, generally, integrable; however, one can assert that it 
also possess N - soliton solutions at any a2 ^ — (l/2)ai, 
which is a generalization of the Hirota solution at q;2 = 0. 
Below, we call Eq. (66) complex MKdV equation. The 
Painleve test [21], applied to Eq. (66), shows that it has 
the Painleve property if 0:2 = 3a3, i.e. 



6a3 = ai = 2a2 



(70) 



This is just the integrability condition for Eq.(l) found 

by Sasa and Satsuma [5] . These conclusions arc in agree- 
ment with those following from the Painleve test of Eq.(l) 
[28,23]. 

If, in addition to condition (3), a2 = (Hirota case 
[4]), we arrive at the complex MKdV equation 

dti^ + 6aM^d^^P + a3dli; = 0. (71) 

The soliton solution to this equation is 

ips{x, t) = a sech[o(a: — ct)\exp[i{px — at)\ (72) 

with arbitrary p and 

c = — 3q:3P^ + 0:30^ , a = —asp^ + 3a3pa^ . (73) 



Then 

*^(X, T) = a sech[o(X - VeT)]ex.p[i{KX - coT)] , (74) 
where 



K+p, Lo = a + n + pV .{75) 



Therefore k is arbitrary and Vs{k), lu{k) coincide with 
Eqs. (2c, d) respectively. This means that Eqs. (74) and 
(75) indeed describe the Hirota solitons. mentioned in 
Sec. 2a. It is also seen that though the Hirota conditions 
do not ensure the integrability ( Painleve tests are neg- 
ative !), Eq. (1) in Hirota case definitely has exact N 
-soluton solutions. 

Let us now consider numerical solutions of the initial 
value problem to Eq.(l) in Hirota case 



ai = 60:3 = —0.6 , Q!2 = 0. 
As initial condition, wc choose a complex pulse 
'1 



(76) 



*(X,0) =^sech 



:(^-^o) 



exp[iC{X - Xo)] (77) 



with A=l , C^O.5 and Xq = 85.3. The structures of 
*(X,r) at T=0 and T=168 are shown in Figs. 4(a,b). 
We see that ^{X,T) evolves into two pulses, well sep- 
arated at T=168 . The spetrum of ^{X,T) at T=168 
(Fig. 5) gives an evidence that, in addition to these pulses, 
no significant wavetrains arc produced at initial condition 
(77). The time variation of the amplitude of the pulse 
with largest amplitude, l^*^^'!, is shown in Fig. 6. It ap- 
proaches the constant limit, approximately equal to 1.5. 
The velocity of that pulse approaches the limit approx- 
imately equal to 0.35 . These results suggest that the 
leading pulse asymptotically approaches a regular soliton 
with a(i) w 1.5 and vi^^ « 0.35. The shapes of «'(X,T), 
obtained numerically for T > 150, confirm this conclu- 
sion. Then, using Eqs. (2c, d) and (4), we can calculate 
the shifted frequency uj' = uj{k) - kvP, which is the 
soliton frequency in the reference frame where the soli- 
ton is at rest. The plot of w' versus V^, shown in Fig. 7, 
gives cj'(0.35) « 1.5. This must be equal to arg^''^^^ On 
the other hand, the numerical solution permits to obtain 
arg^f^^) versus time in the frame where the pulse is at 
rest. This also gives arg^'(^) « 1.5, at sufficiently large 
T . Similar results were obtained for the second pulse. 
Thus we conclude that the initial complex pulse (77) de- 
cays into two Hirota solitons. Qualitatively similar re- 
sults were obtained for the initial pulse 



FIG. 4. 'I' versus a; at ai = -0.6, 02 = 0, aa = -0.1 
and initial condition (77); solid line: |*|, dashed line: Tie*!/ 
(a)T=0. (b) T=168. 
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*(X, 0) = exp 



16 



X-Xo' 



Xo 



85.3 . 

(78) 



FIG. 8. Development of the initial pulse (77) with A=1.9, 
C=0, Xo = 100 (dashed line), (a) ai = -0.5, = 0, as = 
t=48, (b) ai = -0.5, a2 = 0, as = -0.1 t=48, (c) ai = -0.5, 
a2 = 0, as = -0.1 t=175. 

FIG. 5. logu)\'^K\ versus k at ai = —0.6, a2 = 0, 
as = -0.1 and T=168. 



FIG. 6. Temporal behavior of the largest soliton amplitude 
at ai — —0.6, a-z = 0, as = —0.1 and initial condition (77). 

Other features of the dynamics of Hirota solitons, in 
particular their colUsions , will be examined in another 
paper. 



FIG. 9. Plots of Xo{t) for three largest solitons shown in 
Fig.8. Dashed line: ai = —0.5, a2 = 0, as = 0; solid line: 
ai = —0.5, 02 = 0, as = —0.1 



VI. PULSE EVOLUTION AT ai 7^ Gag, = 

Here, we report on the numerical solutions of Eq.(l) in 
general case. First we take, as initial condition, the pulse 
(77) with A=1.9, C=0 and Xq = 100 and assume that 

ai = -0.5 , as = 0.1 . (79) 

The behavior of solution at these parameters is shown in 
Fig.8, where the case 0:3 = (with ai = —0.5) is also 
presented for comparison. In the latter case, the initial 
pulse splits into four solitons of the form (33), (35)- (38), 
with a2 = 1 and moving to the left (Fig. 8a). On the 
other hand, in the case (79) the initial pulse splits into 
three radiating solitons. Two of them propagate to the 
right and the smallest one - to the left. 



FIG. 7. ui' = u! — kVs versus Vs for the largest soliton. 



FIG. 10. Fourier spectra of the solutions shown in Fig.8; 
(a) ai = —0.5, a2 = 0, as = 0; t=48 (b) ai = —0.5, a2 = 0, 
as = -0.1; t=175. 

The radiating wavetrains propagate to the right with 
the group velocities larger than the soliton velocities 
(Figs. 8b and 8c); this is in agreement with Eq. (60). 
The time dependence of the soliton coordinates is shown 
in Fig. 9 . At as, the soliton velocities are constant. At 
aa = —0.1, the solitons are accelerating in positive direc- 
tion, which agrees with Eq. (64). 

The Fourier spectra for both cases, as = and 
as = —0.1, are shown in Fig. 10. The spectrum at 
as = exhibits no resonant radiation (Fig. 10a) while at 
as = —0.1 we see three distinct resonant peaks with neg- 
ative K, which are the wave numbers of the resonant ra- 
diation emitted by the solitons in the "laboratory" frame 
(Fig. 10b). The difference between the wave numbers fol- 
lows from Eq. (20) which is valid, as we have seen, both 
for embedded and radiating solitons. As far as the soli- 
tons have different velocities and amplitudes (the latters 
are determined by the soliton parameter A), they have 
different wavenumbers because k, being a root of Eq. 
(20), depends on the soliton velocity and amplitude. 
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FIG. 11. The amplitude of the first soliton versus time, (a) 
as = 0, (b) as = -0.1. 



FIG. 12. The amplitude of the second soliton versus time, 
(a) as = 0, (b) as = -0.1. 

An important difference between the two cases, as = 
and as ^ 0, is seen in Figs. 11 and 12. In the first case, 
the average soliton amplitudes have constant limits at 
f — > oo while in the second case they must slowly decrease 
because of the soliton radiation. 

In the aforegoing, wc studied the case when a\az > 0. 
Now let us consider a\a-i < 0. We take 



ai = —0.6 , as = 0.1 



(80) 



0:3 > 0, > 0, t/ < 0). On the right hand side, there is 
a radiation composed from the harmonics of continuous 
spectrum; it satisfies the linearized Eq.(l) and has been 
emitted in the transient period of time. The spectal dis- 
tribution is shown in Fig. 14. One can see two spectral 
maximums at positive wavenumbcrs. The narrow one 
corresponds to the resonant radiation, while the broader 
peak is composed of the continuous spectrum. Its struc- 
ture can be understood from the dispersion equation 



uj = (l/2)fc2 - ask^ 



(81a) 



and the correspondig expression for the group velocity of 
continuous spectrum 



(816) 



Prom this we see that, at as > 0, the continues radiation 
on the right hand side of the pulse is composed from 



< A; < fco = l/6a3 



(81c) 



where Vg{ko) = maxyg(fc). As far as ko is less then the 
wavenumber of resonant radiation, approximately given 
by Eq.(50) (note that in our case , according to (45), 
02 > 1), one can see why the peak of continuous spectra is 
less than k of resonant radiation. All this shows that the 
pulse in Fig. 13 is nothing than RRS. The time behavior 
of its amplitude is shown in Fig. 15 ; qualitatively, it is 
similar to Fig. 3a. 

It is reasonable to compare these results with the case 
ai = —0.6, as = , at the same initial condition. Then 
we have two regular solitons, propagating to the right. 
The soliton positions versus time for both cases are shown 
in Fig. 16. Note a very small acceleration of RRS in the 
direction of group velocity (which is now negative) . 



FIG. 13. The RRS at parameters ai = -0.6, as = 0.1. 
Full line: |*|, dotted line : 7?e*. 



FIG. 15. Soliton amplitude versus time at 
ai = —0.6, as = 0.1. 



FIG. 14. Spectral distribution at ai = —0.6, as = 0.1. 

and initial condition (77) with C — 1/2, Xq = 850. 
The solution at T = 370 is shown in Fig. 13. Wc sec a 
pulse with radiation on both its sides. An analysis shows 
that on the left hand side there is the resonantly gener- 
ated radiation (in agreement with Eq. (60), where now 



VII. CONCLUSIONS 

We have considered soliton and quasisoliton solutions 
of Eq. (1) and their relationship. Solutions (2), describ- 
ing regular solitons, degenerate at as 0; this shows 
that they exist due a balance between nonlinear terms 
and linear third order dispersion. On the other hand the 
quasisolitons (embedded and resonantly radiating) 
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FIG. 16. Soli- 
ton positions versus time for qi = —0.6,03 = (1,2) and, 
ai = -0.6, as = 0.1 (3). 

turn at as ^ into the regular sohton solutions of Eq. 
(1) without the third derivative term. Apart from that, 
at a2 = 0, the regular solitons exist only in the special 
case ai = 603. 

The resonantly radiating solitons (RRS) are non- 
steady; the amplitudes of their radiation at small are 
exponentially small and so their parameters change log- 
arithmically slow in this case. Only such quasistcady 
solitons have sufficiently large lifetimes to be of phys- 
ical significance. The embedded solitons (ES), which 
are steady structures consisting of pulses embedded into 
plane waves (wings), have close connection with RRS, 
similar to other systems. The corresponding analytical 
treatment in Sec. Ill is suplemented by numerical simu- 
lations with the cut off operation (Sec. IV) ; it is demon- 
strated that it leads to the transformation of ES into 
RRS that radiates in the direction of the radiation group 
velocity . Just after the cutt off, the amplitude, velocity 
and wavenumber of the radiated wavetrain coincide with 
those of wings; marching in time, we have seen that the 
amplitude decreases, the soliton velocity increases and 
the wavenumber changes according to Eqs.(30) and (31). 
This is in agreement with conservation laws (Appendix 
B). 

In Sec. V we have investigated a special case ai = 60:3, 
interesting from a theoretical point of view . Then the 
quasisolitons, ES and RRS, do not exist and Eq. (1) 
can be reduced to the complex MKdV equation which 
has TV -soliton solutions (consisting of regular solitons). 
This generalizes the Hirota solutions obtained for ai = 
6a3, a2 = 0. The complex MKdV equation possesses the 
Painleve property only at a2 = Sas, which is the case 
when Eq.(l) is integrable. 

Presumably, in the nonlinear processes only regular 
solitons and RRS may be of the physical significance. 
From Eq. (2b) it follows, however, that the amplitude 
of regular soliton vanishes at 0, while for the RRS 

it remains finite. Therefore one should expect that at 
small as, the RRS are more important in the nonlin- 
ear evolution. This was confirmed by solving the initial 
value problem (Sec. VI), where the RRS and not regular 
solitons are seen emerging from the initial disturbances. 
(Small 03 is the most interesting case from physical point 
of view because the third derivative term is in fact the 
result of an expansion) . 



Appendix A 

Modulational instability of plane wave, accord- 
ing to Eq.(l) 

Eq. (1) has exact plain wave solution 

* = Aexp{iKX - iujT) . {AI) 
Substituting (AI) into (1), we have 

w = - (1 - aiK)A^ - a^K^ . {A2) 

Substituting a slightly perturbed wave (AI) 

* = A(l - x)exp(mX - iuiT) , (A3) 
into Eq.(l), we have 

iXT + ^(1 - ^oizK)xxx + i{K + aiA^ - 'iazK'^)xx 

+(1 - aiK)(x + X*)A^ + ia^ixx + Xx)^' • (^4) 

Writing x = u + iw and assuming that 

(u, w) w exp(ipX — irT) . 

we obtain dispersion equation, which is convenient to 
write in the form 

- 2a2kA^T - J(l - 6a3K)k'^[{l - 6a3K)k'^ 

-4(1 - aiK)A'^] = , {A5) 

where 

r = r -I- asp^ - (k - 3a3K^ + aiA^)p . {A6) 

Thus, at real p , Im7 = Imr and the stability condition 
is 

(l-6a3«;)V >4:A^[{l-aiK){l-6a3K)-alA'^] . {A7) 

Unlike the plain wave solution of the regular NLS, now 
the plain wave can be stable at any p , if 

(1 - aiK)(l - 6q!3k) < alA'^ . {A8) 

At a2 = this is possible at 

— <K,<-J— (6a3 < ai) {Ada) 
ai 6a3 

-J— <K< — {ai< 603) {A9b) 
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Applying the stability criterion to the wings of embed- 
ded solitons, where k w 1/2q!3, we see that this may 
satisfy condition (A9b) ; then one can expect that the 
wings are stable. We should also take into account only 
p > L~^, where L is the period in the numerical scheme; 
this relaxes the limitatations following from the stability 
criterion. 

Appendix B 

Investigation of the soliton evolution by means 
of the integrals of motion. 

First , we present a general analysis of conserving in- 
tegrals A^, P, H. Substituting (5) into expressions (12), 
(14), (15) we have after simple algebra 

POO 

N= \^{x,t)\'^dx , {Bl) 

J OO 

1 

P = KN +-i {ipd^i^* - ilj*d^ip)dx . {B2) 

2 Joo 

where 'ip{x,t) satisfies Eq.(9). 

To have a convenient expression for H, we substitute 
in Eq. (15) as^^* and as^^** from Eq.(l), to obtain 

H=- dX{^*dT-^ - ^^T**) + - / dX\^'^ 

^ J OO ^ J OO 

+ -ai / dX\^'^{<i*dx^ - ^dx^*) . 

Then using the Galilei transformation in the form (24) 
and taking into account that at large t 

dt ~ i(AV2)V , 
we arrive at the following asymptotic expression 

Hpi(fl-KV-^X^^N + PV+^qJ \tp\'^dx 
1 f°° 

--iai / m\i^d,r - rd.^)dx . (SS) 

^ Joo 

Substituting (25) and (32) into (Bl), we have at large t 
(cf. Refs. [25,26]), 



N^Ns + Nr, 



where 



/•OO j*00 

Ns= \us{x,t)\'^dx , Nr= \r]{x,t)\'^dx {B5) 

J OO J OO 



are contributions from the soliton and the radiation. Us- 
ing (35) we have 



and, at large t, 



ai 



Nr. 



a2q 



\.f\'\U\t 



{B6) 



m 



where / is given by (52) or (55). Evidently, |/p does not 
depend on x. Then from the conservation of N it follows 



di 



'2a2 ^ ( 




arctan 









-\fm 



{B8) 



In a similar way, from (B2) and the conservation of P we 
have the following asymptotic equation at large t 

dK d f°° 

^Itt^Jtj <i^)d.i^e{x)dx + k\f\''\U\^Q. {B9) 



Prom (36) at a2 = 0, we have 

d^tjjsix) 

Taking into account that 

/ - 

Joo ai \ q 

we transform Eq. (B9) to the form 
^dK d 



"1 2/ \ 



ui^^jdx = ^ ( ^ - AT, 



(SIO) 



{Bll) 



N- 



dt dt 



q_ / 2^/a^A 
ai V q 



A, 



+k\f\'\U\^Q. {B12) 



Now we turn to Eq.(B3). After differentiation over t 
and simple calculations we obtain 



dt 



, d\ ^^dK\ 



d_ 
dt 



-q / \^l^\^dx + -at / |Vr(V'*a,V' - i}d^r)dx 

^ Joo * Joo 



/•OO 
J OO 



{BA) Using 



we obtain 



dV _ dK 
'dt-^'^ltt 



a2(P - KN)— X—N : 

^ ' dt dt 



(S13) 



(S14) 
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1 

dt V 2* 



u^dx + 



2a2 



(B15) 



Let us apply these relations to small ai (which is just 
the case shown in Figs. 3). Assuming that 



aiA 



029 



«1, 



we have 



Therefore Eq. (B8) is reduced to 
"2,/^A afA^ 



d_ 
dt 



1 



= -i/rit^i 



(B16) 



(B17) 



(B18) 



In the same approximation, Eq.(12) takes the form 



^rdK , , ^,2,,,, 2ai d f \^ 



(B19) 



At ai ^ 0, Eqs. (B18) and (B19) turn into equations 
obtained in Rcf.[18]. From (B19), (B14) and Eq. (60) it 
follows an important qualitative result (64). 

Now, consider Eq. (B15). Taking into account 
(B2),(B10), (Bll) and (B19) we see that at 



t <^ xk-^\f\-^u- 



{B20) 



and small ai, Eq. (B15) is approximately equivalent to 
(B18). Performing calculations similar to those in Ref. 
[19] and using (B14), we conclude that the soliton width 
A""'^ and velocity V{t) are changing logarithmically slow. 
The same holds for the soliton amplitude U(9t), deter- 
mined from (39) and (37), in agreement with numerical 
solution of Eq.(l) (see Figs. 3). The chracteristic time of 
the variations of soliton parameters has the same order 
of magnitude as in Ref. [26], which is much less than the 
right hand side of (B20). 

If ai is not small, the analysis is rather tedious and 
will not be considered here. 
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